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ABSTRACT 

We model the overall shape of an accretion disc in a semi-detached binary system in 
which mass is transfered on to a spinning black hole the spin axis of which is misaligned 
with the orbital rotation axis. We assume the disc is in a steady state. Its outer regions 
are subject to differential precession caused by tidal torques of the companion star. 
These tend to align the outer parts of the disc with the orbital plane. Its inner regions 
are subject to differential precession caused by the Lense-Thirring effect. These tend to 
align the inner parts of the disc with the spin of the black hole. We give full numerical 
solutions for the shape of the disc for some particular disc parameters. We then show 
how an analytic approximation to these solutions can be obtained for the case when 
the disc surface density varies as a power law with radius. These analytic solutions for 
the shape of the disc are reasonably accurate even for large misalignments and can 
be simply applied for general disc parameters. They are particularly useful when the 
numerical solutions would be slow. 
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1 INTRODUCTION 

We consider a semi-detached binary system with an 
accretion disc around a black hole the spin axis of 
which is misaligned with the orbital plane of a com- 
panion star. This can arise in a binary system when 
one star has an asymmetric supernova explosion which 
gives a ve l ocity kick to the co ll apsing compact object 
fehklovskiil I l97d : ISutantvol 1 19781 ; iDewev fc Cordesl 1 19871 ; 
iBailes! I198S ). Even a small asymmetry in the super- 
nova can give a large velo city kick to the newly formed 
black hole or neutron star ([Brandt fc Podsiadlowskil 1 19951 ; 
iMartin. Tout fc Pringlejl2009D . 

The tidal torque of the companion causes each ring in 
the disc to precess about an axis parallel to the binary or- 
bital axis. The precession rate increases with increasing ra- 
dius from the black hole. This torque can be balanced by 
the viscosity in the disc to create a steady state disc if the 
orbital period of the companion is short compared with the 
viscous timescale in the outer parts of the disc. The disc 
becomes warped and twisted when the inner regions of the 
disc remain aligned with the spin of the central black hole. 

Microquasars are binary systems in which material is 
accreted from a norm al star on to a compact ob ject with 
relativistic radio jets (|Mirabel fc Rodriguezl 0*9991 ) . The in- 
ner accretion disc aligns with t he spinning black hole by 
the iBardeen fc Pettersonl (1 19751 ) effect. Jets produced in 
this region are expected to be aligned with the spin of 



the black hole. We can measure the inclination of the jets 
and the inclination of the binary orbit and find the mis- 
alignment angle. GRO J 1655-40 might be m i salign ed by 
an ang le of at least 15°. iHiellming fc Rupenl (|l995h mea- 
sured a jet inclination of 85° ± 2° to the line-of-sight and 
iGreene. Bailvn fc Orosj (|200ll ) measured the orbital axis in- 
clination to be 70.2° ± 1.9°. For the mic roquasar V4641 Sgr 
(SAX J1819. 3-2525) . lOrosz et"afl (|200ll ) deduced an orbital 
inclination angle in the range 60° — 70° along the line-of-sight 
and a jet inclination of less than 10° so this microquasar is 
thou ght to be misaligned by at least 5 0°. 

IMartin. Tout fc Pringlel i|200Sf ) and 

IMartin. Reis fc Pringlel (|2008t) examined the microquasars 
GRO J1655- 40 and V4641 Sgr. They used t he steady state 
disc model of IMartin Pringle fc Tout! {2007) warped by the 
IBardeen fc Pettersonl (1 19751 ) effect but included no binary 
torque term so their discs must effectively extend to infinity 
in order to flatten at the outside. Such a disc is aligned with 
the spinning black hole at the centre and at the outside 
the inclination tends to that of the incoming material. We 
consider here what effect a binary torque has on such a 
steady state disc. 



2 WARPED DISC MODEL 

We model the evolution and structure of warped accretion 
discs according to the simplified formalism described by 
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IPringld i|l992t ). We consider the disc to be made up of annuli 
of width dR and mass 2irY,RdR at radius R from the central 
star of mass Mi with surface density ~S(R, t) at time t. The 
angular momentum at radius R is L = (GMi_R) 1/ ' 2 EZ = LI 
where I — (l x ,ly, h) is a unit vector describing the direction 
of the angular momentum o f a disc annulu s with / = 1. 

We use equation (2.8) of lPringld (|l992T ) setting dL/dt = 
and add a term to describe the Lense-Thirring precession 
(if the central object is a black hole) and a binary tidal 
torque term to give 



dL _ 1 d 
~dt ~R~dR 



3R d , rx 3 \ r 1 „ r dl 
-8R^ L) --2^) L+ 2^ RL 8R 



dl 



dR 



L + T LT + Ttid 



The Lense-Thirring precession term is given by 



R 3 



where 



2GJ 



(1) 



(2) 



(3) 



(Kumar fc Pringld [l985h . The angular momentum of the 
black hole J — Jj, with j = (jx,j y ,jz) and \j\ — 1, can 
be expressed in terms of the dimensionless spin parameter 
a such that 



J = ac Mi 



(GMi\ 



(4) 



This approxim ation of the torque is val id if J7lt afi where 
O. 2 = GM/R 3 (|Kumar fc Pringldll985h . This corresponds to 
R 2> R g = GM/c 2 . The inner edge of an accretion disc is 
the innermost stable orbit around the black hole. For a non- 
rotating black hole this is at 6 Jig and decreases down to i? g 
for a maximally rotating black hole. The binary torque term 
is given by 



Ttid = i"Jtid x L. 



(•») 



We discuss this torque further in the next Section. 

There are two viscosities. First vi corresponds to the 
azimuthal shear (the viscosity normally associated with ac- 
cretion discs) and secondly V2 corresponds to the vertical 
shear in the disc which smoothes out the twist. The sec- 
ond viscosity acts when the disc is non-planar. With the 
a-parameterisation the viscosities are given by vi = otiC B H 
and V2 = ct2C s H where c a is the sound speed and H is 
the scale height of the disc. The dimensionless parame- 
ters ai ^ 1 and a.2 must be determined experimentally. 
iPapaloizou fc Pringld (|l983h find that the viscosities can be 
related by 



Vl 



a > 

Oil 



1 

2oJ 



(6) 



to first order. 

The m ore detailed a nalysis of warped discs by lOgilviel 

j 19991 ) and lQgilvid {2000) supports this approach. They in- 
clude, into equation |T}, an extra precession term of the 
form 



i d ( UT . di 



(7) 



where vz is a third viscosity parameter rel ated to the preces- 
sional effects with vz = a^CaH. However. lLodato fc Pringld 
l)2007l ) find that 03 <<iC 02 and we may neglect this term. 

In order to obtain analytic results we take both viscosi- 
ties to have a power law form so that 
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and V2 
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(8) 



where vio, 1^20 and /3 are all constants and Ro is some fixed 
radius. In a steady flat (unwarped) disc the surface density 
far from the inner edge is then 



So 



(9) 



where x — for a steady accretion disc. For a steady decre- 
tion disc in which matter is added at a steady rate at the 
inner radius x = 1/2. This value of x gives the radial depen- 
dence of the surface density at late times in the evolution 
of a decretion disc of fixed mass. An accretion disc reaches 
steady state on the viscous timescale given by 



EL 

Vl 



(10) 



As long as the mass transfer rate remains constant over this 
timescale then the disc can reach this steady state solution. 



3 TIDAL TORQUE 

iLubow fc Qgilvid (|2000h and lOgilvie fc Dubud l|200ll) find 
the tidal torque in the frame of the binary orbit, of radius 
Rh, to be approximately 



T tid 



GM2RT, 
2R 2 



(i) 



3/2 



(e z .l)(e z x I) (11) 



to first order in the angle of tilt between the disc and the 
orbital plane. Here e z is a unit vector in the direction of the 
binary orbital axis, M2 is the mass of the binary companion 
star and i?b is the binary separation. The Laplace coefficient 
of celestial mechanics can be approximated for small 2 by 



6 (1) 



3/2 ( 2 ) = 2 2 
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(12) 



We note that the tidal torque satisfies l.Tud = so that 
the tidal torque is zero when the disc and orbit are aligned. 
Thus even for large tilt angles the tidal torque term is likely 
to have this form, although with modified coefficients. This 
tidal torque is averaged over the orbital period of the sys- 
tem which we assume to be much shorter than the viscous 
timescale at the outer edge of the disc. 



4 RELEVANT DISC RADII 

We expect the centre of the disc to be aligned with the spin 
of the central black hole by the Bardeen-Petterson effect. We 
expect the tidal torques to align the outer regions with the 
binary orbit by an analogous physical process. We note that 
if the outer parts of the disc are closer to counter-alignment 
than alignment, then the whole disc is counter-rotating. The 
inner parts of the disc are flat and counter-aligned with the 
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black hole spin. The counter aligne d disc solution is oppo- 
sitely twisted to an aligned so lution (|Scheuer fc Feilerll 19961 ; 
iMartin. Pringle fc Toutil2007l ). 

The torque on the black hole is always towards align- 
ment of the black hole with the outer parts of the disc be- 
cause the angular momentum of t he disc is effectiv ely very 
large compared to that of the hole (|King et al.ll2005h . This is 
because the angular momentum of the binary orbit is trans- 
ferre d to the outer parts of the di sc when the tidal torque 
acts (|Martin, Tout fc Pringle]|200Sft . 

By comparing the sizes of various terms in the governing 
equation ([TJ we are able to get an initial idea of the structure 
that the disc is likely to have when it is in a steady state. In 
either case, aligned or counter aligned discs, there are three 
radii of physical significance. 



4.1 The Lense-Thirring radius Rut 

A Lense-Thirring radius Rut where the effects of Lense- 
Thirring precess i on an d vis cosity balance was defined by 
IScheuer fc Feilerl <|l996l ) and lMartin. Pringle fc Tout] (|2007V ) 



RtT = 



2UUT 

via Rq' : 



l/(l+/3) 



(13) 



by balancing the first term of equation (TJJ involving vi with 
the Lense-Thirring term and replacing d/dR with 1/R. Well 
inside this radius we expect the disc to be flat and aligned 
with the spin of the central black hole. 



4.2 The tidal radius R ti d 

We define in a similar way a tidal radius, iitid, to be where 
the tidal term balances the same viscous term in the disc so 
that 



.Rtid = 



2v2 Ro P {GM 1 )2Ri 
3GM 2 



2/(7-2/3) 



(14) 



If a disc extends beyond this radius then the viscous ef- 
fects are dominated by the binary torque term. Thus we 
expect that well outside this radius the binary torque term, 
coupled with the viscosity, flattens the disc and aligns with 
the binary orbital plane. The precise form of R tl d will be- 
come clearer in Section f6 . 2 1 when we use it to simplify equa- 
tion ([26)1. 



4.3 The warp radius R wa r P 

The overall warping of the disc is caused by the misalign- 
ment between the Lense-Thirring precession which domi- 
nates at small radii and the tidal precession which domi- 
nates at large radii. The rate of warping of the disc with 
radius is likely to be greatest when these two precessional 
effects balance. We call this the warp radius. It is 



Rw 



3GiM 2 



2/9 



(15) 



It is independent of the magnitude of the viscosity because 
we are now balancing just the two external torques. 



We may then write 



Ri 



+/3 p(7-2/3)/2 
LT -"-tid 



2/9 



(16) 



Thus over all we expect the disc to be warped at radii be- 
tween the Lense-Thirring radius Rut and the tidal radius 
Rtid, with the rate of warping greatest at around the warp 
radius R W arp- 



5 NUMERICAL MODELS 

We are unable to solve the full non-linear disc evolution 
equation analytically and thus in general it is necessary 
to solve the equation s numer i cally. We use the numerical 
scheme described by IPringlei <|l992h . including both a bi- 
nary torque precession term and a Lense-Thirring precession 
term. The binary torque is implemented in the form 



fitid = 



2 (Rud 



R 

Rtid 



3/2 



(17) 



with 7"„ 2 = R 2 /u2, corresponding to just the first term in 
the expansion in equation (|12|l . W e have assumed that th e 
warping is small so that e z .l w 1 (|Qgilvie fc Dubudl200ll ). 
The Lense-Thirring term is 



1 f-g.v 

t^CRlt) \RutJ 



(18) 



We use 200 logarithmically distributed grid points between 
the inner radius Rin and the outer radius R out = 10 4 i?i n . 

We choose two different set of boundary conditions. In 
the first we have zero torque at the boundaries so that E = 
at R = i?i n and R = R ou t which allows mass to flow freely 
through both boundaries. The surface density is chosen ini- 
tially to correspond to a steady state for a flat accretion disc 
with the corresponding accretion rate M , 



E = ( 




if R < i?add, 



Rin 

~R 



if R > R a dd 



(19) 



<|Pringldll992T ). We add mass at i? a dd = 9000 R[ n at a rate 
that keeps the surface density constant at that point with 
angular momentum corresponding to the local Keplerian 
values and aligned with e z . We expect this to be in the 
outer flat part of the disc. Inside this radii the disc acts as a 
standard accretion disc. We set dl/dR = at R = Ri n and 
R — Rout- 

These boundary conditions are realistic but to compare 
directly with our analytic models (Section[6} we use a second 
set of boundary conditions. We compute numerical models 
in which initially E = M/(3ttvi) and the surface density is 
fixed at Ri n and R ou t- The disc is aligned with the black 
hole at the inside so at R — Ri n , I — (sin#o, 0, cos So) and at 
the outside R — R ou t, I = (0,0, 1). 

We choose the viscosities ui — 0.1 (R/Ri n ) 2 and v-z = 
1.0 (R/Rin) 2 so that (3 = 2 and the local viscous timescale 
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r„ = R 2 jv is independent of radius. This minimises compu- 
tational run times. We have chosen v\jvi — 0.1 which cor- 
responds to Qi = 0.2 2 with equation (l6| . With SPH simula- 
tions of warped discs. lLodato fc Pringld (|2007f ) find that the 
ratio V\jv2 is small. Because J7lt oc R~ i and fitid oc R 3 ^ 2 we 
may fix the various relevant radii by choosing the strengths 
of the tidal and Lense-Thirring terms. Models shown in the 
next sections are at a time of 500 (Rf n /GMi ) 2 = 50 = 
500 r„ 2 so that they have reached steady state. 

5.1 Misalignment of 6 = 20° 

We first choose parameters so the Lense-Thirring radius is 
Rut = 10 Rin and the tidal radius is i? t id = 100 -Ri n - We have 
10 R g < Rut < 60 R g and so we can use the approximation 
to the Bardeen-Petterson effect in equation ([3|. For = 2 
the warp radius is then -R warp = 21.5 Ri n (equation \I6\1 . The 
outer disc and tidal precession is aligned with the z-axis. 
The Lense-Thirring precession axis is tilted at an angle of 
80 — 20° to the z-axis. The tilt angle of the disc as a function 
of radius is shown in the left panel of Fig. \T\ The relevant 
radii are also marked. Note that the disc itself is not only 
tilted but that the tilt varies with radius so that the overall 
shape of the disc is a twisted warp (right panel). 

We see that the maximum rate of change of tilt does 
occur at about -R warp . Once we are well within Rut or well 
outside Rtid the disc does indeed become locally flat and 
aligned with the locally dominant precession axis. We see 
that both the numerical models with different boundary con- 
ditions and the analytical model are very similar for this 
small misalignment angle. In Fig. [2] we plot the surface den- 
sity for the two numerical models in Fig. [1] The models are 
similar except in the regions close to the boundaries as ex- 
pected. 

5.2 Changing the Warp Radii 

Here we demonstrate the effect of changing the outer two 
warp radii. We have three models each with the constant 
inclination and surface density boundary conditions. We 
choose Rut and i? t id so that -R warp = 21.5 remains un- 
changed. The first model has the same warp radii as that 
in Fig. Q] with R hT = 10 R in and R tid = 100 R ia . The sec- 
ond has Rut = 3.16 -Rin and Rtid = 10 3 Ri n and the third 
Rut = -Rtid = -Rwarp = 21.5 iiin- We use E = M/(3irvi) at 
-Rin and -Rout- 

In Fig. [3] we plot the inclinations and twists of these 
models. The steepest rate of change of tilt angle occurs at 
around -R wa r P in all models. As .Rtid increases the outer parts 
of the disc become aligned with the outer precession vector 
(orbital plane) at larger radii. As Rut decreases the inner 
disc aligns with the black hole spin at smaller radii. So -R warp 
tells us where the disc is warped Rut and Rtid tell us how 
sharp the warp is. 

5.3 Misalignment of 9o = 60° 

In Fig. [4] we plot the inclination and twist of a disc which is 
misaligned by 60°. The numerical model with the fixed in- 
clination and surface density boundary conditions is similar 
to the analytic solution even at this high inclination. 




12 3 4 



log(R/R bt J 

Figure 2. The surface density for numerical models with black 
hole spin misaligned to the binary orbital axis by 80 = 20°. An 
accretion disc with both the Lense-Thirring effect and a binary 
torque term with Rut = 10 -Rin, Rtid = 10 2 Ri n and R ou t = 
10 4 i?i n - All models have u\, u% oc R 2 . The solid line has the zero 
torque boundary conditions S = at R ln and Rout- The dotted 
line has £ = M"/(37r^i) at R; n and -Rout as in our analytic models. 



5.4 Misalignment of 6 = 80° 

In Fig. [5] we plot the inclination and twist of a disc which 
is misaligned by Go — 80° . We use only the boundary condi- 
tions E = M/(3m>i) at iiin and -Rout- At this high inclina- 
tion the non-linear terms have a bigger effect and the model 
with the same boundary conditions as the analytic solution 
is not such a good fit. 



6 ANALYTICAL MODELS 

To obtain a full solution to the warped disc problem it is 
necessary to resort to numerical means. However it is pos- 
sible to obtain a good analytical approximation to the disc 
structure for 9 ^ 60° . Here we construct such a solution and 
then investigate to what extent it provides a satisfactory fit 
to the full solution. 

The full equation governing the evolution of disc tilt is 
non-linear but, as noted by Scheuer & Feiler (1996), when 
the disc tilt is small enough the non-linear term may be 
neglected. Once the equation has been reduced to a linear 
one, analytical techniques may easily be employed (see also 
Martin, Pringle & Tout 2007). 

For small tilts to the z-axis the local unit vector I — 
(l x ,l y ,lz) that describes the direction of the angular mo- 
mentum can be written to first order as I = (l x , l y , 1) with 
lx, ly <S 1- We then work in terms of the complex variable 

W = l X +ily. (20) 
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Figure 1. Models with the black hole spin misaligned to the binary orbital axis by 9q = 20°. Left: The inclination of the disc. Right: The 
twist of the disc. The models are of an accretion disc with both the Lense-Thirring effect and a binary torque term with Rlt = 10Ri n , 
fttid = 10 2 Ri n and -Rout = 10 4 Ri n . The solid lines have the boundary conditions £ = at R; n and -Rout- The dotted lines are a 
numerical model with £ = M/(3ttvi) at R; n and R ou t which corresponds to the analytic boundary conditions. The dashed line is the 
matched analytical solution described in Section 16.31 All models have v\ , U2 oc R 2 . 




Figure 3. Numerical models of an accretion disc with both the Lcnsc-Thirring effect and a binary torque term with the constant 
inclination and surface density boundary conditions. The black hole spin is misaligned to the binary orbital axis by 20° and vi and 
V2 °c R 2 . All three models have R W arp = 21.5 Ri n - The solid lines have Rlt = 3.16 Ri n and Rtid = 10 3 Ri n - The dotted lines are the 
numerical models in Fig. [TJwith Rlt = 10Rj n and -R t id = 100 Ri n . The dashed lines have R W arp = Rlt = Rtid = 21.5 Ri n . We use 
S = M/(37Wi) at Ri n and Rout- 
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Figure 4. Models with the black hole spin misaligned to the binary orbital axis by 0$ = 60° . The models show an accretion disc with 
both the Lense-Thirring effect and a binary torque term with -Rlt = 10 i?i n , Rtid = 10 2 Rin and Rout = 10 4 Ri n . The solid lines have 
zero torque boundary conditions. The dotted lines are numerical models with £ = M/(3ttu^) at R; n and -Rout- This corresponds to the 
analytic boundary conditions. The dashed line is an equivalent analytical matched model. All models have vi,V2 oc R 2 . 




0.2 0.4 0.6 0.E 



Figure 5. Models with the black hole spin misaligned to the binary orbital axis by 8o = 80°. Models show an accretion disc with both 
the Lense-Thirring effect and a binary torque term with Rlt = 10 Rim Rtid = 10 2 Ri n an< i Rout = 10 4 R ln . The dotted lines are a 
numerical model with £ = M/(3nui) at _R; n and R ou t which corresponds to the analytic boundary conditions. The dashed line is an 
equivalent analytical matched model. Both models have v\, vi oc R 2 . 



We take the scalar product of equation {TJ with I and find 

(21) 



RdR 



This has the solution 



(22) 



where C aC c and Cdcc are constants. The first term is the ac- 
cretion disc term for which fiT, — const and the seco nd term 
is th e decretion disc term which has i/iE k R~2 (jPringld 
Il99lf) . 

For an accretion disc we set L — at R = and L = 
(GMiR) 2 E everywhere so Cdec = 0. For a decretion disc we 
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have Cacc = and so more generally we find 



L = (GMii?)5E 



-0-x 



(23) 



where x = for an accretion disc and x — 1/2 for a decretion 
disc. 



6.1 Solution 

Martin et al (2007) and Scheuer & Feiler (1996) have demon- 
strated how to obtain analytic solutions for the linear warp 
equation with a power law density profile and a single added 
precession term. Here we wish to add two precession terms, 
the Lense-Thirring precession at small radii and the tidal 
precession at large radii. We have not been able to find an- 
alytic solutions for this case. However, as we have seen from 
the numerical results, at most radii in the disc one of the 
precession terms dominates. Thus we may obtain an ap- 
proximate solution to the disc by matching two solutions. 
At small radii we use the solution of Martin et al (2007) 
for a disc with Lense-Thirring precession. At large radii we 
obtain a solution for a disc with only tidal precession. We 
then match these solutions at the intermediate radius -R wa r P 
where the precession terms are of comparable but small mag- 
nitude. 



6.2 Accretion Disc with Tidal Precession 

In steady state without the non-linear term and with only 
the tidal torque term equation (fTJ becomes 



I JL 

RdR 



l p T dl 



(24) 



and, substituting our expression for L (equation I23p we find 

(25) 



1 ii r)l 

iu 20 (GM 1 )-R^- x E R^ 



T - 1 9 

tid ~ RdR 



We add the ^-component of equation (|25|) to i — 
the y-component to get W = l x + il y . Then 



T ttd = ^20 (GMr) s Eo-Rg-^ 



dR 



where 

Ttid = (Tud)x + i(T t id)y 
and the linearised tidal torque is 



T times 



(26) 



(27) 



T tid 



3GM 2 XR 2 
4i? b 



I f — ( — 

Rb 



175 
64 



+ • 



iW (28) 



up to 0(W 2 ) with equation (1 1 2 p . Here Rb is the binar y 
separation and we have used l z ~ 1 l|Ogilvie fe Dubu"sll200 ll ) . 

We only take the first term in the power series. This 
introduces an error which is relatively large in the outer 
parts of the disc. However, our outer boundary condition 
means that the torque itself at the outer edge of the disc is 
zero. We expect the disc to remain aligned with the binary 
orbit until R <C -Rb so the actual error remains small. 



With the binary torque warp radius (equation I14|l we 
simplify equation (|26[) governing the shape of the disc to 



iWy 



2-g-x 



(29) 



i d r 3_ x dw 

V dy [ V dy 

where y = R/Rtid- We solve this equation for accretion discs 
with x = to find the shape of each disc. 

For an accretion disc with x — 0, equation (|29|) becomes 

s dW~ 



iWy 3-f> = A. 

dy 



ti- 



dy 



7-2/j 



to find 



We make the substitution z = y~ 

4 \ 2 a 
7 - 2(3 J dz 

We then set W = y~ 1/4 V = z~^V to find 



iW 



dz J 



(30) 



(31) 



z 2fV dV_ v 

dz 2 dz 



+ i 



(7-2/3) 2 ' "\7-2/3 
This is the modified Bessel equation with solution 

4Vi 



= 0. (32) 



V = Al 



2/3 I 



7-2j3~J • ~ "7=3? \7-2p~J ' (33) 

where I and K are modified Bessel functions of the first and 
second kind and A and B are, possibly complex, constants 
of integration. Substituting for our original variables we find 

4aA 



„, -1/4 

Wad = y 

+BK 



Al 



7-2/3 y7- 

4^1, 



2(3 V 

7-2/3 



7-2(3 



y 



(34) 



We expect the outer edge of the disc to align with the 
binary orbit. Thus we expect W — > as R — > Rb in the 
frame of the orbit. The outer boundary condition should be 
dW/dR = at R = iiout- However when -Rout S> -Rtid then 
we can use W = at R = Rout- 



6.3 Matching the solutions 

We expect the solution Wtid to be valid at large disc radii. At 
small disc radii we make use of the solution given by Martin 
et al (2007) for a disc warped by Lense-Thirring precession. 
The general solution in this case was found to be 



Wlt =r~ 



C K 



2sfi 



4(1+9) 



+D I_ 



2(1+13) 



2Vi 



4(l+/3) 



(35) 



ijMartin. Pringle fc Toul l2007h in the frame in which the 
spin of the black hole is in the z-direction r — R/ Rut- Here 
again I and K are modified Bessel functions of the first and 
second kinds and C and D are, possibly complex, constants 
of integration. A disc subject to both Lense-Thirring and 
tidal torques is dominated by this shape when R <g -R wa r P ■ 
In order to match solutions we need first to transform 
Wur to the frame of the binary orbit, so that the axes of both 
solutions are aligned. Note that we still formally require all 
misalignment angles to be small. To transform this solution 
to the frame of the binary we obtain 



VKlt = cosfloK^LTHsinfWl - |Wlt| 2 + ^(Wlt), (36) 
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where 9q is the angle between the spin axis of the black 
hole and the normal to the orbital plane. We note that we 
are introducing non-linear terms by this transformation but 
because it is only a rotation of the axes it does not affect 
the shape. 

We now take these two solutions, the inner solution Wur 
and the outer solution Wtid, and match them at the inter- 



mediate radius _R„ 



the radius at which the two preces- 



sion terms have an equally small effect. In this way we find 
a matched analytical solution valid at all radii throughout 
whole disc. We are working in the frame of the binary. We 
expect the disc to be aligned with the black hole spin and 
so take 



Wur(Rin) = sin 6> . 



(37) 



This gives one relation between the constants of integration 
C and D. At the outer edge we expect the disc to be aligned 
with the binary orbit and so take 



Wtid(i?out 



0. 



(38) 



This gives one relation between the constants of integration 
A and B. 

At R = -Rwarp we require the solution for W to be con- 
tinuous so that 



WLT (Rwarp) = W tid {R 



(39) 



Because W satisfies a second order differential equation we 
also require that the derivative of W be continuous so that 



dWur 



(Rv. 



dW tic 



(Rv, 



(40) 



dR v w *" p/ dR 
These two conditions give two further relations between the 
constants A, B, C and D so we now have enough to solve for 
the complex constants A, B, C and D and find the complete 
matched solution for the disc. The form of the constants is 
complicated and does not add anything particular to our un- 
derstanding so we don't reproduce them here. We do how- 
ever construct many such matched solutions and compare 
them with our numerical models. 



7 ACCURACY OF THE ANALYTIC 
SOLUTIONS 

We plot these analytical solution in the dashed lines for the 
parameters of the numerical models in Fig. [TJ [3] HI and [5] 
alongside the numerical calculations. At an inclination 9o = 
20° the solutions match very well. Even at 9 = 60° they 
match well when the same boundary conditions are used. 
Though deviations are larger at 80° much of the shape is still 
reproduced. We note that the analytical solutions are easily 
found for any /3. The numerical solution takes much longer 
to evaluate when j3 is small and the timestep, controlled by 
the centre of the disc, becomes very small. 



7.1 Error in neglecting the non-linear terms 

We have assumed that the warp is gradual enough that we 
can neglect the non-linear term l.d 2 l/dR 2 = —\dl/dR\ 2 . 
The relative error in neglecting this term is given by 



Si 



1,. r 8 2 l 



I 



5 7im( L ' lL )- 1 



(41) 



0.01 



0.005 




log(R/Rj 

Figure 6. The effect of neglecting the non-linear terms in the 
matched analytical solution for a disc with the spin axis of the 
black hole misaligned with the binary orbital axis by 8q = 20°. 
We choose Rlt = 10 .Rim Rtid = 100 Ri n and vz/ux = 10 as in 
Fig. [2 



We can justify neglecting this term if E\ <C 1. With similar- 
power- law viscosities 



Ei 



R A 



dl 



dR 



(42) 



We plot this for the matched solution for a misalignment of 
8o — 20° in Fig. |S] We see that the error is at a maximum 
of just over lper cent at R = Rlt- In Fig. [7] we plot the 
maximum value of E\ for different misalignment angles be- 
tween the black hole spin and the binary orbit. It begins to 
dominate after 9o = 60° and this accounts when most of the 
error in the model for 9 = 80°. 



8 SHAPE OF THE DISC 

Given a position in cylindrical polar coordinates (R,(j>i) in 
the x — y plane of the binary orbit, we can find the height, z, 
of the disc at that point. In Cartesian coordinates the disc 
can be described by 



x — R sin 9 cos 4>i 
y — R sin 9 sin <f>i 
z = R cos 9 cos(<j>i — 4>), 



(43) 



where 9(R) — cos _1 (Z z ) is the angle of the direction vector 
of the disc, I, at radius R, to the z-axis. The azimuthal 
angle of I is (p = tan -1 (l y /l x ). In Fig. [8] we picture this 
disc solution, with a logarithmic scale in the radial direction 
and Rin mapped to the origin, using the plotting package 
Povray. We can see the warp and twist of the disc for a 
disc with a misalignment of #o = 45° with the warp radii 
Rut = W Rin and R tid = 100 .Rim 

These models can be used to investigate disc properties 
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Figure 7. The maximum error in neglecting the non-linear terms 
for a model with Rlt = 10 -Ri n and i?tid = 100 Ri n and with 
"a/fl = 10 for a range of misalignments 9q between the black hole 
spin and the orbital angular momentum of the binary system. 

and compare with observations. We illustrate this by exam- 
ining the obscuration of the central black hole by the disc 
at various inclinations. In Fig. [9] we plot the visibility of the 
central black hole on the sky for do = 20° and 60° for the nu- 
merical model with zero torque boundary conditions and for 
the analytically matched solution. The binary orbital axis is 
inclined at i to the line of sight. The angle <j> is the phase 
around the dis c. It changes on the precessio n timescale of 
the black hole dMartin. Tout fc Pringie]|2008l ). This is long 
compared with the orbital period so it is essentially fixed 
relative to the observer. As expected, the higher the mis- 
alignment, the more obscured the black hole is. In each plot 
the region between the two lines would be obscured. Known 
systems would therefore tend to be observationally selected 
to be aligned. 



9 CONCLUSIONS 

We have made both numerical and analytical models of ac- 
cretion discs warped by the Lense-Thirring effect and a bi- 
nary companion. The models have v oc B? with (3 = 2 so 
that the viscous timescale is constant throughout the disc. 
This means that the models converge at their fastest rate. 
More realistically f) = 3/4 and the numerical disc evolves 
slowly because the small viscous timescale at the inner edge 
dominates. On the other hand our analytical matched mod- 
els can be constructed just as quickly for any /3. Up to a 
misalignment of about 9 = 60° they provide a very good 
fit to the disc structure and so may be used where many 
calculations are required in a short time. This is necessary 
when attempting a fit to observations with unknown mis- 
alignment. With a relatively simple analytical fit it is easy 
to compute all sorts of disc properties such as the extent 




Figure 8. A picture of the shape of a disc with the matched 
analytical disc solutions with -Rlt = 10 R; n , Rtid = 100 -R; n and 
a misalignment angle of 8q = 45° . The radial scale is logarithmic. 



of irradiation by the central engine at any point or the in- 
tegrated spectrum without the need to compute a numeri- 
cal model for every set of parameters. We have found that 
known systems are observationally selected to be aligned. 
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Figure 9. The visibility of the black hole with a disc with the numerical (solid lines) and matched analytical disc solutions (dashed 
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to the line of sight and <j> is the phase around the disc. The black hole itself and the central parts of the disc are obscured between the 
two lines (the shaded region). 
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